NASA POWER satellite meteorological system is a good tool for obtaining estimates of the temperature-humidity index under Brazilian conditions compared to INMET weather stations data.

Heat stress negatively affects livestock, with undesirable effects on animals' production and reproduction. Temperature and humidity index (THI) is a climatic variable used worldwide to study the effect of heat stress on farm animals. Temperature and humidity data can be obtained in Brazil through the National Institute of Meteorology (INMET), but complete data may not be available due to temporary failures on weather stations. An alternative to obtaining meteorological data is the National Aeronautics and Space Administration Prediction of Worldwide Energy Resources (NASA POWER) satellite-based weather system. We aimed to compare THI estimates obtained from INMET weather stations and NASA POWER meteorological information sources using Pearson correlation and linear regression. After quality check, data from 489 INMET weather stations were used. The hourly, average daily and maximum daily THI were evaluated. We found greater correlations and better regression evaluation metrics when average daily THI values were considered, followed by maximum daily THI, and hourly THI. NASA POWER satellite-based weather system is a suitable tool for obtaining the average and maximum THI values using information collected from Brazil, showing high correlations with THI estimates from INMET and good regression evaluation metrics, and can assist studies that aim to analyze the impact of heat stress on livestock production in Brazil, providing additional data to complement the existing information available in the INMET database.

Phenotypic variation in milk fatty acid composition and its association with stearoyl-CoA desaturase 1 (SCD1) gene polymorphisms in Gir cows.

Individual variation in milk fatty acid (FA) composition has been partially attributed to stearoyl-CoA desaturase 1 (SCD1) gene polymorphisms in taurine breeds, but much less is known for Zebu breeds. This study investigated the phenotypic variation in milk FA composition, and the influence of SCD1 variants on this trait and on milk fat desaturase indices (DI) in Gir cows. The functional impact of SCD1 variants was predicted using bioinformatics tools. Milk and blood samples were collected from 312 cows distributed in 10 herds from five states of Brazil. SCD1 variants were identified through target sequencing, and milk FA composition was determined by gas chromatography. Phenotypic variation in milk FA composition fell within the range reported for taurine breeds, with SCD18 index showing the lowest variation among the DI. Fourteen SCD1 variants were identified, six of which not previously described. Regarding the A293V polymorphism, all cows were homozygous for the C allele (coding for alanine), whereas all genotypes were detected for the second SNP affecting the 293 codon (G > A), with compelling evidence for functional effects. Significant associations (based on raw p-values) were found between this SNP and C12:0, cis-9, trans-11 CLA and short-chain FA, and between another SNP (rs523411937) and C15:0 and odd-chain linear FA. A new SNP on Chr26:21277069 was associated with trans-11 C18:1, cis-9, trans-11 CLA, C18:3 n-3 and n-3 FA. These findings indicate that SCD1 polymorphisms also contributes to the phenotypic variation in milk FA composition of Gir cows, with potential use in their breeding programmes.

Genomic prediction in Brazilian Guzerá cattle: application of a single-step approach to productive and reproductive traits.

This study aimed to investigate the feasibility of genomic prediction for productive and reproductive traits in Guzerá cattle using single-step genomic best linear unbiased prediction (ssGBLUP). Evaluations included the 305-day cumulative yields (first lactation, in kg) of milk, lactose, protein, fat, and total solids; adjusted body weight (kg) at the ages of 450, 365, and 210 days; and age at first calving (in days), from a database containing 197,283 measurements from Guzerá males and females born between 1954 and 2018. The pedigree included 433,823 animals spanning up to 14 overlapping generations. A total of 1618 animals were genotyped. The analyses were performed using ssGBLUP and traditional BLUP methods. Predictive ability and bias were accessed using cross-validation: predictive ability was similar between the methods and ranged from 0.27 to 0.47 for the genomic-based model and from 0.30 to 0.45 for the pedigree-based model; the bias was also similar between the methods, ranging from 0.88 to 1.35 in the genomic-based model and from 0.96 to 1.41 in the pedigree-based model. The individual accuracies of breeding values were evidently increased in the genomic evaluation, with values ranging from 0.41 to 0.56 in the genomic-based model and from 0.26 to 0.54 in the pedigree-based model. Even based on a small number of genotyped animals and a small database for some traits, the results suggest that ssGBLUP is feasible and may be applied to national genetic evaluation of the breed to increase the accuracy of breeding values without greatly impacting predictive ability and bias.

NAToRA, a relatedness-pruning method to minimize the loss of dataset size in genetic and omics analyses.

Genetic and omics analyses frequently require independent observations, which is not guaranteed in real datasets. When relatedness cannot be accounted for, solutions involve removing related individuals (or observations) and, consequently, a reduction of available data. We developed a network-based relatedness-pruning method that minimizes dataset reduction while removing unwanted relationships in a dataset. It uses node degree centrality metric to identify highly connected nodes (or individuals) and implements heuristics that approximate the minimal reduction of a dataset to allow its application to complex datasets. When compared with two other popular population genetics methodologies (PLINK and KING), NAToRA shows the best combination of removing all relatives while keeping the largest possible number of individuals in all datasets tested and also, with similar effects on the allele frequency spectrum and Principal Component Analysis than PLINK and KING. NAToRA is freely available, both as a standalone tool that can be easily incorporated as part of a pipeline, and as a graphical web tool that allows visualization of the relatedness networks. NAToRA also accepts a variety of relationship metrics as input, which facilitates its use. We also release a genealogies simulator software used for different tests performed in this study.

Genetic study of quantitative traits supports the use of Guzerá as dual-purpose cattle.

OBJECTIVE: The aim of this study was to estimate genetic parameters for 305-day cumulative milk yield and components, growth, and reproductive traits in Guzerá cattle. METHODS: The evaluated traits were 305-day first-lactation cumulative yields (kg) of milk (MY305), fat (FY305), protein (PY305), lactose (LY305), and total solids (SY305); age at first calving (AFC) in days; adjusted scrotal perimeter (cm) at the ages of 365 (SP365) and 450 (SP450) days; and adjusted body weight (kg) at the ages of 210 (W210), 365 (W365), and 450 (W450) days. The (co)variance components were estimated using the restricted maximum likelihood method for single-trait, bi-trait and tri-trait analyses. Contemporary groups and additive genetic effects were included in the general mixed model. Maternal genetic and permanent environmental effects were also included for W210. RESULTS: The direct heritability estimates ranged from 0.16 (W210) to 0.32 (MY305). The maternal heritability estimate for W210 was 0.03. Genetic correlation estimates among milk production traits and growth traits ranged from 0.92 to 0.99 and from 0.92 to 0.99, respectively. For milk production and growth traits, the genetic correlations ranged from 0.33 to 0.56. The genetic correlations among AFC and all other traits were negative (-0.43 to -0.27). Scrotal perimeter traits and body weights showed genetic correlations ranging from 0.41 to 0.46, and scrotal perimeter and milk production traits showed genetic correlations ranging from 0.11 to 0.30. The phenotypic correlations were similar in direction (same sign) and lower than the corresponding genetic correlations. CONCLUSION: These results suggest the viability and potential of joint selection for dairy and beef traits in Guzerá cattle, taking into account reproductive traits.

Identification and in silico characterization of structural and functional impacts of genetic variants in milk protein genes in the Zebu breeds Guzerat and Gyr.

Whole genome sequencing of bovine breeds has allowed identification of genetic variants in milk protein genes. However, functional repercussion of such variants at a molecular level has seldom been investigated. Here, the results of a multistep Bioinformatic analysis for functional characterization of recently identified genetic variants in Brazilian Gyr and Guzerat breeds is described, including predicted effects on the following: (i) evolutionary conserved nucleotide positions/regions; (ii) protein function, stability, and interactions; (iii) splicing, branching, and miRNA binding sites; (iv) promoters and transcription factor binding sites; and (v) collocation with QTL. Seventy-one genetic variants were identified in the caseins (CSN1S1, CSN2, CSN1S2, and CSN3), LALBA, LGB, and LTF genes. Eleven potentially regulatory variants and two missense mutations were identified. LALBA Ile60Val was predicted to affect protein stability and flexibility, by reducing the number the disulfide bonds established. LTF Thr546Asn is predicted to generate steric clashes, which could mildly affect iron coordination. In addition, LALBA Ile60Val and LTF Thr546Asn affect exonic splicing enhancers and silencers. Consequently, both mutations have the potential of affecting immune response at individual level, not only in the mammary gland. Although laborious, this multistep procedure for classifying variants allowed the identification of potentially functional variants for milk protein genes.

Genetic parameters for somatic cell count (SCC) and milk production traits of Guzerá cows using data normalized by different procedures.

This study aimed to estimate the genetic parameters for somatic cell count (SCC) and the genetic association between SCC and milk production traits using two different methods of SCC normalization. The dataset contained information on 8870 lactation records of 6172 Guzerá dairy cows selected for dual-purpose from 95 herds. The lactation means of SCC were normalized in two ways: (a) SCC1 = log10 (SCC) and (b) SCC2 = log2 (SCC/100) + 3. Multivariate analyses were performed considering milk production traits over the course of 305 days of lactation. Estimates of the variance components and genetic parameters were carried out by the Bayesian inference method, applying Gibbs sampling. Single chains of 2,000,000 iterations were used, with sampling discards of the first 5000 chains and a sampling period of every 50 iterations. The deviation of information criteria (DIC) was used to evaluate the best transformation for standardization of the SCC data, comparing analysis 1 (milk production traits over 305 days and SCC1) with analysis 2 (milk production traits over 305 days and SCC2). According to the data structure of this study, SCC1 normalization was the most efficient method and produced better estimates than normalization by the SCC2 method. The heritability estimates for SCC were low regardless of the transformation method used, indicating a small possibility of expressive genetic gains from the direct selection of these traits. However, the repeatability indicated the potential for increasing heritability estimates if the effects of the permanent environment were reduced. The genetic correlations between the milk yield and SCC traits do not indicate the possibility of a correlated genetic gain from the direct selection of one trait. However, concomitant selection for milk production traits and SCC will likely not affect the individual response either.

Genetic parameters for milk, growth, and reproductive traits in Guzerá cattle under tropical conditions.

This study aimed to estimate the genetic parameters of milk (305-day milk yield (MY305)), growth (weaning weight (WW), yearling weight (YW), and weight at 550 days (W550)), and reproductive (age at first calving (AFC)) traits in Guzerá cattle by using Bayesian multiple-trait models. Systematic effects included sex and age at calving for the growth and milk traits, respectively. The additive genetic and contemporary groups (herd and year and season of birth) were included as random effects. Additionally, maternal genetic and permanent effects were also included as random effects for the WW. The heritability estimates were 0.29 (MY305), 0.42 (WW), 0.49 (YW), 0.56 (W550), and 0.25 (AFC). The genetic correlations among the growth traits were higher than 0.83; between the MY305 and WW, MY305 and YW, and MY305 and W550, the genetic correlations were 0.25, 0.32, and 0.36, respectively. The AFC was negatively correlated with the milk and growth traits. These results suggest the viability and potential of the joint selection for milk, beef, and reproductive traits in Guzerá cattle.

Identification of Candidate Genes for Reactivity in Guzerat (Bos indicus) Cattle: A Genome-Wide Association Study.

Temperament is fundamental to animal production due to its direct influence on the animal-herdsman relationship. When compared to calm animals, the aggressive, anxious or fearful ones exhibit less weight gain, lower reproductive efficiency, decreased milk production and higher herd maintenance costs, all of which contribute to reduced profits. However, temperament is a trait that is complex and difficult to assess. Recently, a new quantitative system, REATEST®, for assessing reactivity, a phenotype of temperament, was developed. Herein, we describe the results of a Genome-wide association study for reactivity, assessed using REATEST® with a sample of 754 females from five dual-purpose (milk and meat production) Guzerat (Bos indicus) herds. Genotyping was performed using a 50k SNP chip and a two-step mixed model approach (Grammar-Gamma) with a one-by-one marker regression was used to identify QTLs. QTLs for reactivity were identified on chromosomes BTA1, BTA5, BTA14, and BTA25. Five intronic and two intergenic markers were significantly associated with reactivity. POU1F1, DRD3, VWA3A, ZBTB20, EPHA6, SNRPF and NTN4 were identified as candidate genes. Previous QTL reports for temperament traits, covering areas surrounding the SNPs/genes identified here, further corroborate these associations. The seven genes identified in the present study explain 20.5% of reactivity variance and give a better understanding of temperament biology.

Cluster analysis of breeding values for milk yield and lactation persistency in Guzerá cattle / Análises de agrupamento dos valores genéticos para produção de leite e persistência da lactação em bovinos da raça Guzerá

ABSTRACT: The aim of this study was to explore the pattern of genetic lactation curves of Guzerá cattle using cluster analysis. Test-day milk yields of 5,274 first-lactation Guzerá cows were recorded in a progeny test. A total of 34,193 monthly records were analyzed with a random regression animal model using Legendre polynomials to fit additive genetic and permanent environmental random effects and mean trends. Hierarchical and non-hierarchical cluster analyses were performed based on the EBVs for monthly test-day milk yield, peak yield, lactation persistency, and partial cumulative and 305-day yields. The heritability estimates for test-day milk yields ranged from 0.24 to 0.52. Cluster analysis identified animals in the population that belong to different groups according to milk production level and lactation persistency.

RESUMO: Objetivou-se neste estudo explorar o padrão das curvas de lactação genéticas de bovinos da raça Guzerá, empregando análises de agrupamento. Os 34.193 registros mensais de produção de leite foram provenientes de 5.274 vacas da raça Guzerá, participantes do teste de progênie. As análises foram realizadas com um modelo de regressão aleatória com polinômios de Legendre, composto pelos efeitos aleatórios genético aditivo, de ambiente permanente e o residual, e a curva média de lactação da população. Análise de agrupamento hierárquico e não hierárquico foram realizados com base nos VG para a produção acumulada até os 305 dias, pico e persistência da lactação, e períodos parciais da lactação. As estimativas de herdabilidade para produção de leite no dia do controle variaram entre 0,24 a 0,52. A análise de agrupamento identificou os animais da população que pertencem a diferentes grupos de acordo com o nível de produção de leite e persistência da lactação.

Novel SNPs and INDEL polymorphisms in the 3'UTR of DGAT1 gene: in silico analyses and a possible association.

Diacylglycerol-O-acyltransferase (DGAT1) gene encodes the rate-limiting enzyme of triglyceride synthesis. A polymorphism in this gene, DGAT1 K232A, has been associated with milk production and composition in taurine breeds. However, this polymorphism is not a good tool for ascertaining the effects of this QTL in Bos indicus (Zebu), since the frequency of the DGAT1 232A allele is too low in these breeds. We sequenced the 3'-untranslated region of DGAT1 gene in a sample of bulls of the breeds Guzerá (Bos indicus) and Holstein (Bos taurus) and, using in silico analysis, we searched for genetic variation, evolutionary conservation, regulatory elements, and possible substitution effects. Six single nucleotide (SNPs) and one insertion-deletion (INDEL) polymorphisms were found in the Guzerá bulls. Additionally, we developed a preliminary association study, using this INDEL polymorphism as a genetic marker. A significant association was detected (P ≤ 0.05) between the INDEL (DGAT1 3'UTR INDEL) and the breeding values (BV) for protein, fat, and milk yields over a 305-day lactation period. The DGAT1 3' UTR INDEL genotype I/I (I, for insertion) was associated with lower BVs (-38.77 kg for milk, -1.86 kg for fat, and -1.48 kg for protein yields), when compared to the genotype I/D (D, for deletion). I/D genotype was lower D/D genotype (-34.98 kg milk, -1.73 kg fat, and -1.09 kg protein yields). This study reports the first polymorphism of DGAT1 3'UTR in the Guzerá breed, as well as its association with BV for milk protein, fat, and milk yields.